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INTRODUCTION 


Many  important  physical  systems*  exhibit  highly  unstable  or  chaotic  behavior.  Traditional  methods 
for  treating  such  systems  have  been  linearization  and  stability  analysis.  It  is  generally  acknowledged  that 
non-linear  ^sterns  exhibit  important  features  that  cannot  be  understood,  controlled,  or  exploited  without 
addressing  ^eir  behavior  in  a  non-linear  form.  Thus,  techniques  are  being  developed  to  model  and 
analyze  sudi  complex  systems.  The  approach  entails  the  anafysis  of  complicated  structures  that  exhibit 
fractal  geometry.  Fractal  parameters  describe  how  the  length,  area,  and/or  density  of  an  oi^'ect  varies  with 
scale  and  nuy  provide  quantitative  measures  of  the  complex  structures  that  occur  in  nature.  To  optimally 
exploit  fractal  perspectives,  quantitative  techniques  for  measurement  of  fractal  properties  are  evolving.  In 
this  report,  qualitative  and  quantitative  techniques  are  applied  to  elucidate  the  fractal  structure  of  the 
strange  attractors  that  govern  the  evolution  of  symmetric  chaos  in  some  model  systems. 

BACKGROUND 

The  qualitative  features  of  the  time  evolution  of  many  physical  systems  are  embodied  in  the 
solutions  to  equations  of  the  form: 

t)  =  F(x(t) ,  a,  b,  . . . )  (la) 

where  a,  b,  etc.  are  parameters  and  x  e  R”.  The  function  F  determines,  in  principle,  the  non-linear 
evolution  of  the  system.  In  practice,  the  behavior  of  non-linear  dynamical  systems  is  quite  complex,  and 
the  solutions  to  Eq.  (la)  (i.e.,  the  trajectories  x(t))  are  frequently  extremely  sensitive  to  small  changes  in 
initial  conditions  or  the  system  parameters.  Equation  (la)  is  a  continuous  time  description  of  the 
dynamical  system,  however,  such  systems  are  general^  treated  in  terms  of  difference  equations  or  discrete 
time  mappings  of  the  form 

x{/i  +  l)  »  F(x(n)  ,  x(n-l)  ,  . . . ,  a,  b,  . . . )  (li?) 

The  motion  of  non-linear  systems,  governed  by  Eq.  (lb),  may  be  unbounded  or  may  settle  down 
after  an  initial  "transient"  to  a  subset  of  R”  called  an  attractor.  The  motion  of  such  systems  can  be 
characterized  by  their  attractors  rather  than  by  specific  trajectories.  The  attractors  of  non-linear  dynamical 
systems  fail  into  two  categories: 

i.  Periodic  attractors.  When  an  attractor  comprises  a  finite  set  of  points  (that  are  visited  in 
sequence),  the  motion  of  the  system  is  periodic  and  one  refers  to  a  periodic  attractor.  If  the  number  of 
points  is  N,  it  is  called  a  period  N  attractor. 

ii.  Chaotic  attractors.  When  an  attractor  contains  an  infinite  set  of  points,  the  motion  of  the 
system  is  called  "chaos,"  trajectories  are  called  "chaotic,"  and  one  refers  to  "strange  attractors"  or  "chaotic 
attractors."  Strange  attractors  have  fractal  structure. 

The  attractors  studied  here  describe  the  dynamics  of  an  important  subset  of  chaotic  systems: 
symmetric  chaotic  mappings  in  E^  which  were  qu^tatively  characterized  by  Field  and  Golubitsky^.  In 
particular,  we  focus  our  efforts  on  the  quantitative  determination  of  the  fractal  parameters  that 
characterize  the  attractors  of  mappings,  which  transform  from  D.toZ.  symmetry  as  a  parameter  is 
changed.  Field  and  Golubitsky^  referred  to  the  motion  of  non-linear  systems,  whose  attractors  were 
invariant  with  respect  to  such  symmetry  groups,  as  "symmetric  chaos." 

Box-counting  algorithms  for  fractal  analysis  are  based  on  a  statistical  description  of  the  variation 
with  scale  of  the  distribution  of  "mass"  of  a  fractal  set.  We  employ  a  box-counting  technique  to  extract  the 
fractal  parameters  that  characterize  the  attractors  of  four  systems  exhibiting  symmetric  chaos. 
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APPROACH 


The  attractors  of  the  potynomial  mappings  determined  fay  iteration  of 
z  -  {k  a.  \z\^  *  B  Re{z°)  +i(o)z  +  Y  {z*)°  (2) 

where  z  e  C,  and  A,  a,  B,  /,  and  a  are  real  constants,  have  D,  symmetry  for  a>sO,  and  7^  symmetry  for 
I  a>  I  >0.  iMeld  and  Golubitsky^  recently  discussed  qualitative  techniques  for  anatyzing  the  attractors  of 
potynomial  mappings. 

The  nature  of  the  attractors  of  Eq.  (2)  is  elucidated  by  generating  graphical  representations  of 
their  structure,  similar  to  those  of  Ref.  2,  and  by  fractal  analy^ 

The  purpose  of  the  graphical  representations  is  to  visually  illustrate  the  chaotic  structure  and 
symmetry  of  the  attractors  studied.  Although  a  number  of  subtle  decisions  must  be  made  in  such 
representations,  only  standard  techniques  are  employed. 

The  fractal  analysis  provides  quantitative  measures  of  the  attractor  structures.  The  fractal  analysis 
employs  the  agglomeration  box-counting  (ABQ  algorithm  of  Meisel,  Johnson,  and  Cote^.  ABC  is 
especially  suitable  for  analysis  of  large  point  sets. 

Since  the  ABC  algorithm  is  not  a  standard  procedure,  we  outline  the  technique  here.  ABC 
enqtloys  a  set  of  hypercubes  (boxes)  of  linear  dimension  E  to  cover  a  subset  of  the  fractal  set  For  each  E, 
the  partition  function 

Z(<7;£)  (3) 

1 

is  computed  for  a  range  of  the  descriptive  index  q,  and  p;(E)  is  the  occupation  probability  for  the  i***  box  of 
side  E.  The  partition  function  is  related  to  the  Hentschel  and  Procaccia'^  fractal  dimension  D(q)  by  the 
scaling  relationship: 

Z(q;E)  ~  £-<9- 1)  (4) 

where  D(q)  is  a  generalization  of  the  Mandelbrot  fractal  dimension^  -At  qsQ,  0(q)  reduces  to  the 
standard  Mandelbrot  fractal  dimension.  If  D(q)  is  non-integer  and  invariant,  the  set  is  monofractal  If 
D(q)  varies  with  q,  it  is  said  to  be  multifractaL 

The  essential  steps  in  applying  ABC  are: 

i.  Define  a  set  of  elementary  boxes  for  the  given  point  set  S,  such  that  the  elementary  boxes  have 
edge  lengths. 


E(H,K,M,...)=E<,/(2«  3*'  5“...) 

where  H,  K,  M, ...  are  integers,  and  Eg  is  to  be  the  largest  box  edge  length  in  the  set. 

ii.  Compute  (measure,  etc.)  and  store  the  occupation  numbers  (integers)  n,(E(H,K,M,...))  for  the 
elementary  boxes.  (The  probabilities  p,(E)=n,(E)/N  need  not  be  stored.) 
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iiL  For  each  h  <  {0,1,..^},  k  <  {0,l.»MK},ni  e  define  sets  of  boxes  of  edge  length 

E(h.k^...)»Eo/(2*  3"  5-...), 

which  also  cover  the  point  set  and  contain  integer  numbers  of  elementary  boxes. 

a.  Compute  the  occupation  numbers  nt(E(h,k4n,...))  for  each  box  by 
summing  the  occupation  numbers  in  the  elementary  boxes 
contained. 

b.  For  a  representative  set  of  q  values,  compute  Z(q,£(h,k,m,...))fiom 
Eq.  (3). 

iv.  For  each  q  value,  compute  D(q)  by  least  squares  fitting  of  a  straight  line  to 
ln(Z(q,E(hJMn,...))/(q-l)vs  In(E(hJun,...))  as  follows  from  Eq.  (4).  In  applications,  the  elementary  box 
size  must  be  large  enou^  to  contain  a  "substantial  number”  of  members  of  the  set,  and  the  largest  boxes 
included  in  the  fit  must  be  smaller  than  the  extent  of  the  point  set  (Special  techniques  are  required  for 
qa>l,  but  D(q)  is  a  smooth  function  of  q  in  any  case.) 

ABC  has  been  tested  on  Euclidean  point  sets  and  sinq>le  monofiractal  and  multifractal  constructs 
for  which  anafytic  values  of  D(q)  could  be  computed.  Worst  case  results  were  within  5  percent  of  the 
analytic  values  for  q20  and  10’  for  the  cases  studied  in  Ref.  3.  Reliable  values  of  D(q)  were  not 
obtained  employing  ABC  (or  conventional  box-counting*)  for  q<0;  therefore,  D(q)  for  q<0  are  not 
reported  here. 

RESULTS 

Results  are  presented  for  four  classes  of  symmetric  chaotic  mappings  based  on  Eq.  (2):  Figures  1 
and  2  pertain  to  attractors  having  3-fold  symmetry  (i.e.,  D,  and  2^  synunetry).  Fig.  3  to  attractors  having  6- 
fold  symmetry,  and  Fig.  4  to  attractors  having  16-fold  symmetry. 

For  each  class  of  mapping?  defined  by  a  specific  choice  of  A,  a,  B,  and  n,  attractors  were 
generated  based  on  10’  iterations  of  Eq.  (2)  for  (o  values  starting  from  0  and  stepping  by  0.01  over  the 
range  for  which  the  motion  is  bounded.  Smaller  changes  of  of  were  taken  in  "interesting"  ranges.  The 
parameter  set  defining  each  class  of  mappings  is  given  in  the  captions  of  the  a-parts  of  each  figure.  (The 
A,  a,  B,  y,  and  n  parameter  sets  were  chosen  to  match  the  Field  and  Golubitsky’  mappings.) 

Two-dimensional  images  of  subsets  of  the  attractors  (generated  for  at—O  and  selected  other  values 
of  (w),  which  were  formed  using  the  density  of  the  attractor  mapping  as  an  index  to  a  gray  scale  using  a 
technique  similar  to  Field  and  Golubitsky^,  are  shown  for  each  case.  The  attractors  are  depicted  in  square 
regions  whose  comers  are  at  (-2,  -2)  and  (2,  2)  for  the  3-fold  mappings  and  (-1,-1)  and  (1,1)  for  the  others. 
Each  image  is  based  upon  768  x  768  pixels  shaded  according  to  number  of  hits  modulo  512. 

Three-dimensional  images,  in  which  the  smoothed  density  of  the  attractor  is  the  third  coordinate, 
are  shown  for  the  6-  and  16-fold  symmetric  attractors  to  enhance  visualization.  Althou^  a  more  faithful 
representation  was  obtained  prior  to  smoothing,  unsmoothed  density  is  not  shown  because  the  rapid 
(chaotic)  density  variations  were  difficult  to  interpret. 

Graphs  of  the  Hentschel-Procaccia’  generalized  fractal  dimension  D(q)  (for  q=0,  2,  5,  and  9)  vs 
over  the  range  for  which  chaotic  attractors  exist  and  D(q)  vs  q  for  the  ca  values  corresponding  to  the 
chaotic  attractor  images  are  shown  for  all  cases. 
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For  aQ  the  mappings  presented: 

L  For  a»l,  unbounded  motion  results. 

ii.  For  ai<l, 

a.  a  range  of  values  of  at  for  which  periodic  orbits,  possibly 
comprising  parts  of  period  doubling  cascades  with  decreasing  (u, 
was  found; 

b.  chaotic  attractors  were  found  for  a  range  of  positive  (o  extending 
from  0; 

c.  dramatic  changes  in  the  structure  of  the  attractors  with  varying  a> 
were  apparent  in  their  two*  and  three-dimensional  images. 

D(0)  is  essentially  invariant  under  changes  in  to  except  for  the  mappings  of  Fig.  1;  however,  D(q) 
exhibits  substantial  non-monotonic  variation  with  to  over  the  chaotic  range  for  aU  classes  of  mappings 
studied. 


The  3-fold  symmetric  mappings  of  Figs.  1  and  2  are  most  interesting: 

The  mappings  of  Fig.  1  yield  a  single  basin  of  attraction  at  small  to  that  fissions  into  three 
symmetrically  disposed  sub-basins  of  attraction  for  to  near  033.  Near  the  "critical  point”  one  passes 
through  a  "period  window."  An  enhanced  rendering  (i.e.,  each  point  is  rendered  as  a  small  circle)  of  one 
of  the  three  period  8  sub-attractors  in  the  periodic  window  at  a>a<0338S00...  is  shown.  An  image  of  a 
corre^jonding  chaotic  sub-attractor  at  to»Q39  is  also  shown.  The  full  attractors  in  the  three  sub-basin 
range  have  2^  symmetry.  However,  the  motion  of  the  dynamical  system  would  be  conGned  to  one  of  the 
sub-basins  of  attraction  and  would  not  have  Z,  symmetry.  Searches  based  on  changes  in  cu  of  0.001  did 
not  yield  periodic  orbits  with  a  period  other  than  8  (=2*)  in  the  window.  (Although  the  ABC  algorithm 
yielded  values  as  large  as  0.06,  D(q)  values  are  shown  as  0.00  for  to  in  the  strictfy  periodic  windows  in  Figs. 
1  and  2.) 

The  mappings  of  Fig.  2  yield  a  single  basin  of  attraction  or  unbounded  orbits  for  almost  all  to. 

TVvo  strictly  periodic  windows  were  found,  and  their  attractors  are  shown.  The  window  at  i;vs0.06  contains 
a  period  2  attractor  in  the  midst  of  a  one  basin  chaotic  range  of  to.  The  window  at  a>ai0.12  contains  a 
period  39  (a>3*i3)  attractor,  which  has  2^  symmetry,  but  the  factor  of  13  is  puzzling.  The  window  at 
to»028  contains  a  quasi-period  12  attractor,  which  has  Z,  symmetry,  but  is  not  strictly  periodic  even  after 
10*  iterations.  The  ABC  algorithm  yields  values  of  D(q)  greater  than  030  in  accord  with  the  apparent 
chaotic  nature  of  the  attractor  for  (us038. 

The  mappings  of  Figs.  3  and  4  yield  substantial  variations  in  structure  and  multifractal  D(q)  with 
<».  Although  prominent  dips  were  apparent  in  D(q)  near  specific  values  of  to,  no  periodic  windows  were 
discovered  searching  via  steps  of  0.001  in  to  near  those  dips. 
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CONCLUSIONS 


The  results  presented  in  Figs.  1  through  4  comprise  a  qualitative  and  quantitative  solution  (except 
for  the  possible  omission  of  some  narrow  periodic  windows)  of  the  non-iinear  pr'^blems  defined.  One 
needs  both  approaches  to  characterize  the  dynamics. 

Even  the  simplest  non-iinear  dynamical  systems  exhibit  unexpectedly  conq)lex  behavior.  Indeed, 
even  in  retrospect,  it  seems  impossible  to  deduce  the  complex  form  of  the  symmetric  chaotic  attractors  or 
the  diversity  of  D(q)  vs  q  for  representative  ot-values,  such  as  those  illustrated  in  Figs.  1  through  4,  by 
analysis  of  the  form  of  the  mappings  of  Eq.  (2). 

The  nature  of  the  change  in  symmetry  from  Dg  to  Z,  is  generally  not  reflected  in  the  Mandelbrot 
fractal  dimension  D(0),  but  it  is  manifested  in  changes  of  the  Hentschel  and  Procaccia^’^  generalized 
fractal  dimension  D(q).  Thus,  both  qualitative  and  multifractal  characterization  of  the  attractors  of  non¬ 
linear  systems  of  interest  are  recommended. 
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Rgure  2a.  Attractor  with  n=3;  A=-2.38;  a=1.0;  B=0.0;  y=0.1. 
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Rgure  3b.  Three-dimensional  representation  of  attractor  in  3a 
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Rgure  4a.  Attractor  with  n=16;  A=2.39:  a=-2.5;  B=s-0.1;  y=0.9. 
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Rgure  4b.  Three-dimensionai  representation  of  attractor  in  4a 
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